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Abstract 

A review is given of quantum network models in class C which, on a suitable 
2d lattice, describe the spin quantum Hall plateau transition. On a general class 
of graphs, however, many observables of such models can be mapped to those of 
a classical walk in a random environment, thus relating questions of quantum and 
classical localization. In many cases it is possible to make rigorous statements about 
the latter through the relation to associated percolation problems, in both two and 
three dimensions. 
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1 Introduction 



Lattice models of spatially extended systems have a long record of usefulness in condensed 
matter physics. Even when the microscopic physics is not necessarily related to a crys- 
talline lattice, it can be very useful to concentrate the essential degrees of freedom onto 
a regular lattice whose length scale is larger than the microscopic one yet much smaller 
than that the expected scale of the physical phenomena the model is designed to address. 
In many cases, the phenomenon of universality ensures that this idealization can never- 
theless reproduce certain aspects exactly. The classic example is that of a lattice gas, 
where a coarse-grained lattice on the scale of the particle interaction radius is introduced 
and used to make predictions for continuum systems, in cases in which the correlation 
length is large, for example close to the liquid-gas critical point. 

The lattice models discussed in this article - in this context called network models 
- were first introduced by Chalker and CoddingtonpQ as a theoretical model for non- 
interacting electrons in two dimensions in a strong transverse magnetic field and in the 
presence of disorder: the physical setting for the integer quantum Hall effect. The starting 
point is to consider non-interacting electrons moving in two dimensions in a disordered 
potential V(r) and a strong perpendicular magnetic field B. We assume that the length 
scale of variation of V(r) is much larger than the magnetic length. In this limit the elec- 
tronic motion has two components with widely separated time scales[2]: the cyclotron 
motion and the motion of the guiding center, along contours of V(r). The total energy 
of the electron in this approximation is E = (n + \)fiw c + V(r), where u c = eB/m and 
n labels the Landau levels, and we therefore expect to find extended states at energies 
E c corresponding to those values of V at which the contours of V(r) percolate. In 2d 
this is expected to occur at one particular value of V, which can be taken to be V = 0. 
Otherwise, away from the percolation threshold the guiding centers are confined to the 
neighborhoods of the closed contours, corresponding to bulk insulating phases which con- 
duct only along the edges of the sample. This immediately provides a simple explanation 
of the experimental result that extended states occur only at the transition and not in the 
Hall plateaux. If this were literally correct the plateau transition would be in the same 
universality class as classical percolation. 

However, this picture is modified for energies close to the transition, since quantum 
tunnelling is expected to be important [3] where closed contours approach each other, see 
Fig. [TJ The network model idealizes this picture by distorting the percolating contours 
of V(r) into a regular square lattice, known as the L-lattice, shown in Fig. |5J In this ap- 
proximation the potential V takes a checkerboard form, being > on (say) even squares 
and < on odd squares. In the limit of large magnetic field the spin degree of freedom of 
the electrons can be ignored and there is a one- dimensional vector space associated with 
each edge. The quantum tunnelling at each node is taken into account by an S-matrix 
which connects the spaces on the adjacent incoming and outgoing edges. This depends 
on E in such a way that for E > E c tunnelling between edges bordering regions with 
V > is enhanced, and vica versa for E < E c . Apart from this, the most important 
quantum feature of the problem is the phase which the electron wave function picks up 
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Figure 1: A typical configuration of contours of the random potential V(r). In the 
guiding center approximation the particle follows these. Only those with V ~ (thicker 
lines) are important for the plateau transition. Quantum tunnelling can occur close to 
the saddle points of V(r) (circled). Figure adapted from Ref. pQ. 
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Figure 2: The idealization of Fig.[T]on the L-lattice. The edges correspond to the contours 
with V ~ and tunnelling occurs at the nodes. 

on traversing a closed contour, which is proportional to the magnetic flux through the 
loop, and therefore its area. On idealizing the loops to a regular lattice, this is repre- 
sented by a quenched uncorrelated random flux through each plaquette, or, equivalently, 
a quenched random U(l) phase as the particle traverses a given edge. While in principle 
the connectivity of the lattice, the S'-matrix elements at the nodes, and the U(l) phases 
on the edges are all quenched random variables, in fact only the latter appear to be the 
most relevant in describing the universal properties of the transition. 

The Chalker-Coddington model was initially analyzed numerically [T] using transfer 
matrix methods. Its predictions appear to agree remarkably well with experimental re- 
sults, perhaps embarrassingly so since it ignores electron-electron interactions which may 
become important near the transition. However it has so far resisted all attempts at 
an analytic solution (as have other more sophisticated field theoretic formulations of the 
integer quantum Hall plateau transition [4]), which by now perhaps elevates this to being 
one of the outstanding unsolved problems of mathematical physics. 

Later, following interest in various forms of exotic superconductivity, it was suggested 
that certain disordered spin-singlet superconductors, in which time-reversal symmetry is 
broken for orbital motion but Zeeman splitting is negligible, should exhibit a quantum 
spin Hall effect, in which the role of the electric current is replaced by that of a spin 
current. [5] The single-particle hamiltonians for such a system then turn out to possess 
an Sp(2) (or equivalently SU(2)) symmetry. In the classification scheme of localization 
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universality classes due to Altland and Zirnbauerj6] they are labelled as class C. The 
corresponding variant of the Chalker-Coddington model is straightforward to write down, 
and was studied, once again numerically [TJ [5]. Then, in a remarkable paper, Gruzberg, 
Read and Ludwigjjg] argued that several important ensemble-averaged properties of this 
model (including the conductance) are simply related to those of critical classical 2d per- 
colation. This is a powerful result because many of the universal properties of percolation 
are known rigorously [9j. It therefore gives exact information about a non-trivial quantum 
localization transition. 

The arguments of Gruzberg et al. [8j were based on a transfer matrix formulation of 
the problem and therefore restricted to the particular oriented lattice (the L-lattice) used 
by the Chalker-Coddington model, which is appropriate to the quantum Hall problem 
in 2d. They also used supersymmetry to perform the quenched average. Subsequently, 
Beamond, Cardy and Chalker[10] gave an elementary, albeit long-winded, proof of their 
main result which holds for any lattice of coordination number four, and any orientation 
of this lattice as long each node has two incoming and two outgoing edges. Later, this was 
shown in a slightly more elegant fashion using supersymmetry [TT] . In each case certain 
quenched averages of the quantum problem are related to observables of a certain kind 
of classical random walk on the same lattice. If the quantum states are localized, the 
corresponding classical walks close after a finite number of steps. If the quantum states 
are extended, in the classical problem the walks can escape to infinity. 

Since this correspondence holds on a very general set of graphs and lattices, it can 
be used more generally to improve our understanding of quantum localization problems. 
In particular it can test the generally accepted notion that in two dimensions all states 
are localized, except in certain cases with special symmetries (such as at the Hall plateau 
transition). It can help understand why in higher dimensions there should be in general 
a transition between localized and extended states, and possibly illumine the nature of 
that transition. More mathematically, it may shed light on the search for a rigorous proof 
of the existence of extended states. Apart from the L-lattice considered by Gruzberg et 
al, it is possible in several cases to use known information about percolation to place 
bounds on the behavior of the classical walks and hence the quantum network model. Up 
until recently, these arguments have been restricted to two dimensions, but now suitable 
three-dimensional lattices have been identified in which the correspondence to classical 
percolation is explicit. 

The layout of this paper is as follows. In Sec. [2] we describe general network models 
and observables which are related to experimentally measurable quantities. In Sec. [3] 
we summarize the supersymmetric proof[TT] of the main theorem which relates suitable 
quenched averages of these observables in the Sp(2) network model on a general graph 
to averages in a classical random walk problem on the same graph. The next section 
H] describes how these classical models on certain lattices (the L-lattice, relevant to the 
quantum Hall effect, and the Manhattan lattice) relate to 2d classical percolation which 
can then be used to bound their behavior. In Sec. we extend this to some special 3d 
lattices, describing relatively recent work, some of it so far unpublished. Finally in Sec. 
we discuss some outstanding problems. 
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Figure 3: A 3 — > 3 node which can be decomposed into 2 — > 2 nodes. 

2 General network models 

In this section we define a general network model on a graph Q and discuss the kind of 
observables we would like to calculate. The graph Q consists of nodes n and oriented 
edges. Initially suppose that Q is closed, that is every edge connects two nodes, and that 
each node has exactly two incoming and two outgoing edges. In fact the general theorem 
to be proved in Sec. [3] holds for more general graphs, but it can be shown [TTj that the 
corresponding classical problem has non-negative weights (and so admits a probabilistic 
interpretation) if and only if each node in Q and its correspond transition amplitudes can 
be decomposed into a skeleton graph with only 2 — >• 2 nodes. See Fig. [31 

On each edge e of Q is an N- dimensional Hilbert space H, e . We assume these are all 
isomorphic. The Hilbert space of the whole system is then ® e ^g}i e . We consider a single 
particle whose wave function at at time t is a superposition of the basis states in this 
space. The dynamics is discrete: if the particle is at the center of edge e at time t, at the 
next time t + 1 it must move in the direction of the orientated edges through a node to 
the center of a neighboring edge e'. Because of the discrete dynamics we should consider 
the unitary time-evolution operator U. This has an off-diagonal block elements U e i e which 
have the form 

Ue'e = U^/ 2 S e ' e U^ 2 . 

Here U e is a unitary N x N matrix which maps T-L e — > H e and describes the rotation of 
the wave function in the internal space as the particle moves along the edge e, and S e / e 
maps "H e — > 1-l e > and describes the transmission through a node. The evolution matrix 
after time t therefore has non-zero block elements 

7( e t> e o) 

where the sum is over all Feynman paths j(e t , eo) on Q of length t, starting at and 
ending at e/. Note that in this sum a given edge can be traversed an arbitrary number 
of times. 

Both the matrices U e and S e > e are quenched random variables, assumed chosen from 
the invariant measure on some subgroup of U(iV). This is chosen according to the sym- 
metry class under consideration. For spinless, or spin-polarized, electrons, where electric 
charge is conserved, we can take N = 1 and the U e G U(l). For models with class 
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C symmetry [6], corresponding for example to the spin quantum Hall effect, the single- 
particle Hilbert space is even- dimensional, and there is an action of a y such that the 
single-particle hamiltonian % satisfies %* = ~a y T-La y . This implies a symmetry between 
states with energies ±E, and that the time-evolution operator U l = e~ Mt satisfies 

= OyliOy , 

implying that the matrices U e should be symplectic, in Sp(iV), which for iV = 2 is iso- 
morphic to SU(2). 

For a given node n with incoming edges (ei,e 2 ) and outgoing edges (e'^e^), the S- 
matrix has the block form 

where each element is an N x N matrix. However, since all the matrices are chosen at 
random, there is the gauge freedom of redefining S e > e — > S e ' e V e , U e — > V^UeVe, which 
allows us to choose each S e > e to be proportional to the unit N x N matrix ljy. The 
remaining 2x2 matrix can further be chosen to be real and therefore orthogonal. Thus 
in fact ([T]) can be replaced by 

V — sin U n cos o n I 

We are left with the gauge-transformed U e and the 8 n as quenched random variables. How- 
ever, we shall treat them differently, first keeping the 8 n fixed while performing quenched 
averages over the U e . In fact we shall see that in most cases is suffices to take all the 6 n 
as fixed and equal on each sublattice. 

Because we consider discrete rather than continuous time evolution, the usual Green 
function is replaced by the resolvent (1 — zU)~ l of the unitary evolution operator, whose 
matrix elements we shall however continue to refer to as the Green function: 

G(e f ,e i; z) = (e'\(l-zU)- l \e). (3) 

Here |e) € H, has non-zero components only in T-L e . Note that G(ef,ei]z) is an iV x iV 
matrix mapping 'H ei — > "H e/ . The parameter z is the analog of the energy (roughly 
z ~ e lE ). For \z\ 1, the expansion of ([3]) in powers of z gives G as a sum over Feynman 
paths from ej to e/. Each path 7 is weighted by z' 7 ' times an ordered product of the U e 
with e G 7 and the factors of cos# n or ±sin# n for n G 7. For a finite closed graph Q, this 
expansion is convergent for \z\ < 1 and therefore defines G as an analytic function in this 
region. In general, for a finite Q, there are poles on the circle \z\ — 1 corresponding to 
the eigenvalues of U. 

However, for \z\ > 1, G admits an alternative expansion in powers of z^ 1 by writing 
it as 

G(e f , ei ;z) ee -{e'\z- l U\l - z- l U ] Y l \e) . (4) 

This is given by a sum over paths 7 with length > 1 of a product of 2H 7 ! with ordered 
factors Ul along the path, and cos n or ± sin 9 n as before. 
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The eigenvalues of U have the form e l6j , where the — ir < 6j < it with j = 1, . . . , M 
are discrete for a finite graph Q. We define the density of states by 

j 

In the standard way, the density of states is given by the discontinuity in the trace of the 
Green function, this time across \z\ — 1 rather than ImE = 0: 

p{e) = £ *Z ( Tr G ^ e ' z = ^ ~ Tr G ^ e ' z = e ^)) > ( 5 ) 

where the trace is in the N- dimensional space T-L e . In the case where Q is a regular lattice, 
in the thermodynamic limit we expect the eigenvalues to be continuously distributed 
around the unit circle. 

We note that in the U(l) case, when the U e are pure phases e 1 ^, in each term in the 
Feynman path expansion of fl3Bl) a given edge e occurs with a weight e tne ^ e , where nj 
is the number of times the path 7 traverses this edge. On averaging a given path 7 will 
contribute to the mean density of states only if nj = 0, that is it has length zero. Thus, 
in the U(l) network models, G(e, e, z) = 1 for \z\ < 1, and zero for \z\ > 1, and the mean 
density of states is constant, and completely independent of the 6 n . This is consistent 
with the general result that at the plateau transition in the charge quantum Hall effect, 
the density of states is non-singular. 

We now turn to the conductance. In order to define this, we must consider an open 
graph, which can be obtained from a given closed graph Q by breaking open a subset {e} 
of the edges, relabelling each broken edge e as e in and e out . External contacts are subsets 
Ci n and C out of these. The transmission matrix T is a rectangular matrix with elements 

T=(e out \(l-U)- 1 \e in ), 

where e- lTL G Ci n , e out G C out . Note that for an open graph the resolvent (1 — zlX)~ x 
generally has poles inside the unit circle \z\ — 1. In the thermodynamic limit, however, 
as long as the fraction of broken edges is zero (for example if we have contacts only along 
part of the edge of the sample), it is believed that the limit as \z\ — > 1 can be taken. 
The multi-channel Landauer formula then gives the conductance between the contacts 

as 

g = (Q 2 /h)TiT^T, 

where Q is the quantum of charge carried by the particle. For the integer quantum Hall 
effect, Q = e, and for the spin Hall effect Q = \h. 

3 The main theorems. 

In this section we focus on the Sp(2) (=SU(2)) case and summarize the method of proof of 
the main theorem relating the quenched averages of the density of states and conductances 
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Figure 4: The two ways to decompose a node of Q. Note that in general this does not 
have to be planar, although for the L-lattice (Fig. [2]) it is. 

to observables of a classical random walk model on the same graph Q. We restrict attention 
to graphs with exactly 2 incoming and 2 outgoing edges. The general case is considered 
in Ref. [11] and is considerably more verbose. 

Let us first define the corresponding classical problem. Starting with a given closed 
graph Q, to each node with incoming edges (ei, e 2 ) and outgoing edges (e[, e' 2 ) associate its 
two possible decompositions ((e'^ei), (e' 2 ,e2)) and ((e' 2 ,ei), (e'^e^), corresponding to the 
two distinct ways of passing twice through the node without using a given edge more than 
once, irrespective of the order in which the edges are used. This is illustrated in Fig. HI 
To each decomposition assign a probability p n = cos 2 6 n or 1 — p n = sin 2 9 n , corresponding 
to the ^-matrix in (J2]). Doing this for each node in Q gives a decomposition of Q into a 
union of closed loops. There are 2' e l such decompositions, where \Q\ is the total number of 
nodes, and the assigned probabilities give a product measure on the set of decompositions. 

We now state the two main theorems. We use the symbol A to denote a quenched 
average in the original quantum network model. Let P(L, e) be the probability that the 
edge e lies on a closed loop of length L in the classical problem. 
Theorem \.G{e^e^z) vanishes unless Cf = e^, when it is given by 

rrTT, r_ / 2-£ L>0 P(L,e)z 2L : \z\ < 1 

(e ' e '"° I E L>0 P(L,e)z-^ :\z\>l " 

If we apply this to the density of states in §5§ we find simply 



where P(L) = J2 e e), the probability that a edge chosen at random lies on 

a loop of length L. Note that if this distribution decays sufficiently fast for large L, p(e) 
is expected to be analytic, while if it decays as a power law, p(e) will have a power law 
singularity at e = 0. 

For an open graph with external contacts C- lTL and C out , the decomposition of Q can 
also lead to open paths connecting edges in C- m to those in C out . In this case, for the 
conductance, we have 
Theorem 2. The mean conductance is 



where P(e out , e in ) is the probability that an open path from e in to e out exists. That is, g is 
twice the mean number of open paths connecting C m to C out . 




(6) 



L>0 
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We describe the proof of these two Theorems using the supersymmetric path integral 
method of Ref. [11]: the combinatorial method in Ref. [10] is perhaps more illustrative of 
why the result holds, but it involves heavier algebra. 

In the standard way, the Green function G, being the inverse of a matrix, may be 
written as a gaussian integral over commuting (bosonic) variables. The notation is a little 
cumbersome but the basic idea is simple. Label each end of a given edge e by e R and e^, 
in the direction of propagation e R — > e L . Introduce complex integration variables b R {e) 
and &l(c), each of which is a 2-component column vector in the SU(2) space. Then 

G(e f , er, z) = (b L (e f )b L (erf) = J Ue [db L (e)][db R (e)]e^ (?) 

where W[b) = W edge + W node and 

W edge = z b[(e)U e b R (e) , W node = Y,Y1 b R^i) S sM^i) ■ 

e n ji 

We use the notation (• • • ) to denote averages with respect to this gaussian measure. The 
measure for each integration is 

— 6t6 



[db] = (1/tt 2 ) J e- b]b dRebdlmb. 

For a finite graph there are a finite number of integrations and the integral is convergent 
for \z\ < 1. 

The next step is to average over the quenched random matrices U e . As usual this 
is difficult because these appear in both the numerator and denominator of ([7]). This 
can be addressed using replicas, or, much more effectively in this case, by adding an 
anticommuting (fermionic) copy (/, /) of each pair bosonic variables (b,tf). Note that 
each / is also a 2-component column vector in SU(2) space, and each / a row vector. The 
Grassman integration over these is defined by 



J[df} = J e-^dfdf. 



where 



J [df]f = J [df]f = and j [df]l = J [df]ff = 1 . 
Integrating over these cancels the denominator in (j7J) so that 

G(e f , ei ;z) = [ n[^L(e)][^(e)][rf/ L (e)][^(e)]/L(e / )/ L (e J )e M/ [ b ] +H/ W 



The action W[b] + W[f] is supersymmetric under rotating bosons into fermions and vice 
versa, so we can replace the bosonic fields at and ef by fermionic ones and consider 

{/,.(' f )f,.(< ,))■ 
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Now consider the average over U on a given edge. This has the form 

J dUexp(zb[Ub R + zf L Uf R ) . (9) 

Because the anticommuting fields square to zero, the expansion of the exponential in 
powers of the second term terminates at second order. The gauge symmetry discussed 
earlier shows that the purely bosonic zeroth order term is in fact independent of b\ and 
b R and is in fact unity. The second order term is also straightforward: carefully using 
the anticommuting property we see that is actually proportional to det U — 1, times the 
determinant of the SU(2) matrix of fermion bilinears with elements fufitj, where i and 
j take the values 1 or 2. The first order term then simply converts this into something 
supersymmetric. The conclusion is that the integral ([9} equals 1 + \z 2 det M where the 
2x2 matrix M has components Mij = b^ Li b R j + fufiij- This can be rearranged in the 
form 

1 + z 2 [(l/V2)(b{J L2 - b{J L1 )} [(1/V2)(b m f m - b R2 f m )] + z 2 [f L J L2 ] [f R2 f m ] . 

Each expression in square brackets is an antisymmetric SU(2) singlet. The three terms 
above describe the propagation of either nothing, a fermion-boson (fb) singlet, or a 
fermion-fermion (//) singlet along the edge e. This is a remarkable simplification: before 
the quenched average, each edge may be traversed many times, corresponding to the prop- 
agation of multi-particle states. It is this which gives one of the principle simplifications 
of the SU(2) case, that does not happen for U(l), one of the reasons this is much more 
difficult. 

This result shows that single fermions or bosons cannot propagate alone, so G(e/, e^) = 
if Cf 7^ Ci. A non-zero correlation function is however 

(fi(e f )f2(e f )f 2 (e i )fi(e i )) = G U G 12 - G 21 G 22 = detG(e f ,e u z) (10) 

Now G, being a real linear combination of products of SU(2) matrices, can in fact always 
be written in the form AG where A is real and G GSU(2)H Hence det G = A 2 and 
G^G = A 2 1, so TrG^G = 2 det G. When z — 1 this gives the point conductance between 
ti and tf. 

When Ci = Cf = e, however, we can always insert a pair of fermion fields f 2 {e)f 2 (e) 
into the correlator (/i(e)/i(e)) at no cost, so that in fact 

G(e, e; z) n = G(e, e; z) 22 = det G(e, e; z) . 

Thus both the mean density of states and the mean conductance are proportional to det G 
and therefore are given by the correlation function (jTOl) . 

The next step is to consider propagation through the nodes. Note that we can now 
drop the distinction between b R and bi, etc. The contribution from a given node takes 
the form 

n^( e ^n<^)' ( n ) 

* 3 

2 This follows from the representation U = cos a + i(cr ■ n) sin a. 
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where A 1 = 1, A 2 = (l/>/2)(6a/ 2 - fc/i), A 3 = /i/ 2 , and 

5 = exp ( h \ s vh + fiSijfj 

where SV,- is the matrix in (j2]). Since this expression conserves fermion and boson number, 
it follows that the total numbers of (fb) and (//) singlets are also conserved. Also, only 
terms second order in the Sij survive. In fact, after a little algebra (lllj) reduces td_| 

$ a'^! $ a' 2 a 2 S 11 $22 ~ 5a' 2 a 1 5a' 1 a 2 S2lSl2 = 5a' 1 o 1 ^c 2 cos2 ®n + § o^ot-fi a! x a.% srr l 2 ®n ■ 

These two terms correspond to the decomposition of the node n described earlier. It 
shows that, on performing the quenched average, the quantum network model is equivalent 
to a classical one in which Q is decomposed into disjoint loops, and along each loop 
propagates either an (fb) singlet, an (//) singlet, or nothing. Theorems ( l3|3i) now follow 
straightforwardly. We argued above that the mean diagonal Green function G(e, e; z) is 
given by the (//) correlation function ([/i(e)/2(e)][/i(e)/2(e)]). So, in each decomposition 
of Q, there must be an (//) pair running around the unique loop containing e. This gets 
weighted by a factor z 2L . Around all the other loops we can have either an (//) pair, a 
(bf) pair, or just 1. The (//) pair, being itself bosonic, gives z 2L for a loop of length L, 
while the (bf) pair, being fermionic, gives — z 2L . These two cancel (by supersymmetry) , 
leaving a factor 1 for every loop other than the one containing e. The argument for 
Theorem [3] for the conductance works in the same way. 

Note that these methods may be extended to the quenched averages of other observ- 
ables in the quantum model, although the density of states and the conductance are most 
important. However not all quantities of interest can be treated in this fashion. For 
example, the fluctuations in the conductance involve 



(G(e out ,e in )tG(e out ,e in )) , 

and, in order to treat this, we would need to double the number of degrees of freedom in the 
integral representation. Many of the formulas which are special to SU(2) integrations then 
no longer hold. In this context it is important to note that the conductance fluctuations 
in the quantum model are not given by the fluctuations in the number of paths connecting 
the two contacts in the classical model. (If this were the case, the quantum system would 
be behaving completely classically!) 

An amusing application [10] of the general theorems is to consider single edge e, closed 
on itself, but take U e G Sp(iV) with iV even and > 2 in general. A general Sp(iV) 
matrix may be built up in terms of successive Sp(2) rotations in overlapping 2-dimensional 
subspaces. For example for Sp(4) we may write a general matrix in the form 



_1_ ( U x 0\ / 1 1 

72 V U 2 1 



(12) 



3 In Ref. this was carried out for a general node of arbitrary coordination, with excruciating 
algebra. 



10 



Xx>oo<-XOC 



Figure 5: Graph corresponding to a single link in the Sp(4) model, and its topologically 
distinct decompositions. Each node corresponds to the S"-matrix which is the second 
factor in (fl2l) . so each term in the decomposition is equally weighted. 

where U\ and U2 are independent Sp(2) matrices. If these are drawn from the invariant 
measure on Sp(2), then the product of a large number of independent such matrices will 
converge to the invariant measure on Sp(4). Thus, in a particular Sp(2) basis, Q has 
the form shown in Fig. [51 After applying Theorem [31 each decomposition corresponds to 
a permutation of the different channels corresponding to the basis chosen in f[T2"|) . This 
generalizes to arbitrary N. If we now connect opposite ends of Fig. to make a closed 
graph, we find, after the decomposition, all possible lengths L of loops from 1 to |iV, with 
equal probabilities. Thus P{L) = 2/N for 1 < L < ^N, and zero otherwise. Using (J6]), 
this gives for the density of eigenvalues of a random Sp(iV) matrix 

P{e} ~ 2nN V (iV + l)sin e yl ' 
in agreement with Ref. [T2] . 

4 Two-dimensional models 

In this section we discuss the consequences of the main theorems for specific 2d lattices 
relevant to physical problems. 

4.1 The L-lattice 

This is the lattice, illustrated in Fig. [21 used by the original Chalker-Coddington model for 
the quantum Hall plateau transition. The reasons for choosing this lattice were discussed 
in Sec. [TJ In the class C version of this, the same lattice is used, the only difference being 
that the quenched random U(l) phases on each edge are replaced with SU(2) matrices. 
We also recall that, because of the checkerboard nature of the potential, where even 
plaquettes correspond to V > and odd ones to V < 0, the angles 9 n , which represent 
the degree of anisotropy of the tunnelling at the nodes, are in fact staggered: 9 n = 9 on 
the even sublattice, and (ir/2) — 9 on the odd sublattice. Thus for 9 = all the loops in the 
decomposition of Q will be the minimum size allowed, surrounding the even plaquettes, 
and for 9 = n/2 they will surround the odd plaquettes. Away from these extreme values, 
the loops will be larger. If there is a single transition it must occur at 9 = %/A. 

The mapping to square lattice bond percolation for this model is exact. Consider 
independent bond percolation on the square lattice rotated by 45° with respect to the 
original, whose sites lie at the centers of the even plaquettes of the original lattice. See 
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V 

Figure 6: The mapping between decompositions of the L-lattice and bond percolation on 
the lattice CI . Open bonds are shown as thick lines, closed as dashed lines. 

Fig. [6] Each edge of C intersects a node of the original one. We declare it to be open, with 
probability p = cos 2 9, or closed, with probability 1—p — sin 2 9, according to the way the 
node is decomposed in the classical loop model on Q. There is thus a 1-1 correspondence 
between decompositions of Q and bond percolation configurations on C . We can also 
consider the dual lattice C" whose vertices are at the centers of the odd plaquettes of the 
original lattice. Each edge of this lattice crosses a unique edge of and we declare it to 
be open, with probability 1 — p, if the corresponding edge of CJ is closed, and vice versa. 
For each percolation configuration on £', there is a corresponding dual configuration on 
£," . The clusters and dual clusters fill the plane without intersecting each other. (For this 
purpose it is necessary to regard single isolated sites as clusters.) Close to p c — ~ many 
clusters nest inside dual clusters and vice versa. 

For a given decomposition of Q, the loops give the hulls of percolation clusters on C 
and dual clusters on C" . For a closed simply connected lattice Q, these are closed curves 
which simultaneously circumscribe a cluster and inscribe a dual cluster, or vice versa. For 
open boundary conditions, some of these curves may begin and end on the boundary. 

From this mapping to 2d percolation, and known exact and conjectured results about 
the latter, many results about the SU(2) network model on the L-lattice may be deducedjH]. 
At the transition, the conductance between two bulk points a distance r apart is given 
by the probability that they are on the same loop, which is known to decay as |r|~ 2:ri 
where x\ = ■? [13]. The conductance of a rectangular sample with contacts along oppo- 
site edges is given by the mean number of hulls which cross the sample. At the critical 
point this depends only its aspect ratio Li/L 2 , in a complicated but calculable way|14j. 
For L\fLi ^> 1 (where L\ is the length of the contacts) it goes like g(Li/L 2 ) where g 
is the universal critical conductance which, from conformal field theory results applied 
to percolation, takes the value \/3/2.[H] In Sec. [3] we showed that the mean density of 
states is given by the probability P(L) that a given edge is on a loop of length L. At 
criticality, this probability scales like L~ Xl where df = 2 — x\ = | is the fractal dimension 
of percolation hulls. [13J. Also \z — 1| oc \E\, which is conjugate to L, has RG eigenvalue 
yi = 2 — x\. This means that the singular part of the mean density of states behaves 
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Figure 7: The Manhattan lattice. 




Figure 8: The decomposition of each node of the Manhattan lattice and its relation to 
percolation on CJ . Open bonds shown as solid thick lines, closed as dashed. 

like [8] 

p(E) ~ \E\ Xl/yi = \E\ l/7 . 

A number of other exponents, including the usual percolation correlation length exponent 
v — |, were identified in the physics of the spin quantum Hall transition in Ref. [8]. 

4.2 The Manhattan lattice 

Although the L-lattice is the natural candidate for studying the spin quantum Hall tran- 
sition, the mapping discussed in Sec. [3] is of course valid for any orientation of the edges 
of a square lattice, and one can legitimately ask whether other possibilities lead to in- 
teresting physics. The expectation, based on the continuum classification of localization 
universality classes, is that unless there is some special symmetry, such as occurs for the 
L-lattice with its sublattice symmetry, all states in 2d will be exponentially localized and 
therefore almost all loops in the classical model will have finite length. This was studied 
for the Manhattan lattice in Ref. [15]. On this lattice all edges in the same row or column 
are oriented in the same direction, and these alternate, see Fig. [7[ This resembles the 
one-way system of streets and avenues in Manhattan. At each corner, the driver can go 
straight on, or turn either left or right according to the parity of the intersection. 

Consider an SU(2) network model on this lattice. The decomposition of the lattice 
corresponds to replacing each node either by a crossing, with probability 1 — p, say, or 
an avoidance, with probability p (see Fig. [8]). Note that in this case the loops on the 
decomposed lattice are in general non-planar. Nevertheless it is possible to make rigorous 
progress using a mapping to percolation, owing to the sublattice structure. Consider once 
again bond percolation on the 45°-rotated lattice £ (see Fig. [S]). An edge is declared 
open if the corresponding node is decomposed in such a way that the paths turn by 90°, 
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as shown. The open edges once again form connected clusters which enclose, and are 
enclosed by, percolation hulls, and similarly for the dual clusters. A little reflection shows 
that each loop of the decomposed Manhattan lattice is constrained to lie on or between 
neighboring hulls which enclose a dual cluster. For p > p c = | the dual clusters are almost 
surely finite, and therefore so are their surrounding hulls, and therefore also the loops of 
the decomposed Q. Therefore for p > \ the SU(2) network model on the Manhattan lattice 
is in the localized phase. An unproven conjecture, consistent with our expectations for 
Anderson localization in 2d, is that this happens for all p > 0. Simulations [15] of the 
classical loop model indicate that this is the case for p > 0.2. The field theory arguments 
to be discussed later indicate that there should always be a finite localization length, 
diverging as £ ~ exp (const. /p 3 ) as p — > 0. 

4.3 Other 2d lattices 

We may consider other orientations of the square lattice. In general a given orientation 
corresponds to a configuration of the six- vertex model, which satisfies the 'ice rule' that 
there are two incoming and two outgoing arrows at each vertex, or node. The number of 
such allowed configurations grows exponentially with the size of the lattice, but we could, 
for example, consider a randomly oriented lattice in which the weights for different types 
of node are given by the 6- vertex model. In that case the L-lattice and Manhattan lattice 
are just particular extreme points of the parameter space. Once again we can associate an 
edge of the lattice £ or £' with the decomposition of a node where the path has to turn. 
These edges can be thought of as two-sided mirrors, reflecting all the paths which impinge 
on either side. The study of these 'mirror models' as models for classical localization has 
been extensive (see, for example Ref. [16J, and references quoted therein), although on 
the whole only simulational results are available. However it is important to realize that 
arbitrary mirror models do not in general lead to a unique orientation for the edges: a 
path may traverse a given edge in both directions. Thus the set of mirror configurations 
corresponding to quantum network models is a subset, and it would be wrong to infer 
general conclusions about these from the study of the wider problem. Nevertheless, one 
expects that, for a sufficiently high density of randomly oriented mirrors, the paths are 
finite and so the states are localized. An interesting and unresolved question, however, is 
what happens at low mirror density, when the mean free path is large. Expectations from 
quantum localization would then suggest that the paths in such models are still localized 
on large enough scales, unless there is some special symmetry like that of the L-lattice. 

5 Three-dimensional models 

We now discuss some results for class C network models on 3d lattices. It is possible, 
of course, to consider layered 2d lattices which might be used as models for quantum 
Hall physics in multilayered systems. For example, in a bilayer system consisting of two 
coupled L-lattices, depending on the strength of the coupling between the layers, one 
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expects to see either two separate transitions between states of Hall conductance 0, 1 and 
2. These can be simply understood in terms of the classical model and an equivalent 
percolation problem. |17j However, such models do not capture some of the important 
properties of real bilayer systems which depend on electron-electron interactions. 

One motivation for studying truly 3d class C network models is to shed light on 
the physics of the localization transition in 3d. So far, these have been carried out 
either numerically or by mapping to 3d percolation. However the restriction to four-fold 
coordination in order that the equivalent classical model has non-negative weights [IT] 
means that it is necessary to use loose-packed lattices. 

5.1 Diamond lattice 

The most extensive numerical simulations have been carried out on the diamond lattice. [TS] 
Although this lattice has cubic symmetry, assigning the orientations of the edges breaks 
this down to tetragonal, inducing an anisotropy. However, numerical tests show that this 
is not very great, giving, for example, a ratio of about 1.1 between the conductances in 
the two distinct directions in the conducting phase. As expected, the model exhibits a 
sharp transition at p = p c between an insulating phase and a conducting phase. This is 
shown, for example, in data for the conductance G(p, L) of a cubic sample of linear size L, 
which, according to the Theorem [3j is given by the mean number of open paths between 
two opposite faces. In the insulating phase p < p c this should approach zero as L — > oo, 
while for p > p c we expect ohmic behavior with G(p, L) ~ a(p)L. Close the critical points 
we expect finite-size scaling of the form 

G(p,L) = f(L/t(p)), 

where the localization length ~ \p — p c \~ u ■ Thus the data should show collapse when 
plotted as a function of (p — p^L 1 ^, and this is clearly exhibited in Fig. |9j The best fitted 
value for v, taking into account corrections to scaling, is[18j v = 0.9985 ± 0.0015. The 
closeness of this value to unity, the value predicted by a first-order result u~ l = e + 0(e 2 ) 
of the 2 + e-expansion (see Sec. ?? and Ref. [21]) is remarkable, but perhaps a coincidence. 

At p — Pc, the weighted number of return paths of length L, behaves as P(L) ~ 
L~ xi / yi , as discussed in Sec. 14.11 In Ref. |18j this exponent is denoted by 2 — r, where 
numerically r = 2.184 ± 0.003. This exponent is related to the fractal dimension df = 
y 1 = 3 — X\ of the paths at p c , giving df = 2.534 ± 0.009. By the same arguments as in 
Sec. 14. X ^ this gives, for example, the singular part of the density of states p(E) ~ |£ , | r ~ 2 . 
Note that, because of the mapping to the classical problem for which far larger systems 
can be studied, the error bars on these exponents are much smaller than those quoted for 
the conventional 3d Anderson transition. 

5.2 3d L-lattice and Manhattan lattice 

It is possible to construct 3d oriented lattices with cubic symmetry which are direct 
analogues of the 2d L-lattice and Manhattan lattices discussed in Sec. S] and for which 
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Figure 9: Conductance as a function of (p — p^L 1 ^, illustrating scaling collapse. Upper 
inset: same data on a logarithmic scale. Lower inset: conductance as a function of p for 
several values of L. Lines are a scaling fit described in detail in Ref. [18]. (Reproduced 
from Ref. [18] with permission of the authors.) 

the arguments relating the classical models to percolation can be generalized. 

Consider two interpenetrating cubic lattices C\ = Z 3 and C 2 = (Z + I) 3 . Each face of 
C\ intersects an edge of C 2 at its midpoint. The four faces of C 2 which meet along this edge 
intersect the given face of C\ along two mutually perpendicular lines, also perpendicular 
to the edge (see Fig. [10]). These lines form part of the lattice Q. The same is true, 
interchanging the roles of C\ and C 2 . The full graph Q is the lattice formed by the 
intersection of the faces of C\ with those of C 2 . The nodes of Q lie on the midpoints of the 
edges of C\ (the centers of the faces of C 2 ) and vice versa, and have coordination number 
4. Clearly Q has cubic symmetry. 

For the L-lattice on Q, the orientation of the edges is chosen so that each node looks 
like the nodes of the 2d L-lattice, as in Fig. [101 There is an overall 2-fold degeneracy 
in assigning these, but once the orientation at one node is fixed, so are the rest. The 
S'-matrices of the network model, and the corresponding probabilities p and 1 — p for 
the decomposition of Q are assigned consistent with the percolation mapping now to be 
described. 

The sites of C\ may be assigned to even and odd sublattices C[ and C'{ according to 
whether the sum of the coordinates is even or odd. Each of these lattices is in fact a 
face-centered cubic (fee) lattice. Now consider nearest neighbor bond percolation on C[. 
Each nearest neighbor edge of C[ intersects the midpoint of an edge of C 2 , along which 4 
faces of C 2 intersect. A decomposition of this edge consists in connecting up these faces 
in neighboring pairs, as in Fig. [TTJ 

For each edge there are two possible decompositions, and we can do this in such a way 
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Figure 10: A node of Q (purple) is formed by the intersection of four faces of C 2 (green) 
and one face of C\ (red), and vice versa. The orientation shown corresponds to the 3d 
L-lattice. 




Figure 11: Decomposition of the four faces (in red) at an edge of C\ corresponding to 
bond percolation (bonds shown in brown) on an fee sublattice of C 2 . 



that if the corresponding edge of the percolation problem on C[ is open, it passes between 
the connecting pairs (see Fig. ITTj) . just as in Fig. E]in 2d. If on the other hand the edge 
is closed, then it intersects both pairs. Equivalently, we can consider percolation on the 
'dual' fee lattice C'[. Each edge of this lattice intersects one edge of C[ at the midpoint 
of an edge of C 2 . We declare the edge of C" to be open if the intersecting edge of C[ is 
closed, and vice versa. 

For a finite lattice, each decomposition of the edges of C 2 divides the faces of C 2 into 
a union of dense, non-intersecting, closed surfaces, in the same way that a decomposition 
of the corresponding nodes of the 2d square lattice divides the edges into non-intersecting 
closed loops. These closed surfaces form the hulls of the bond percolation clusters on the 
fee lattice C[ and its dual C". That is, each closed surface either touches a unique cluster 
externally and a unique dual cluster internally, or vice versa. 

However, this is only half the description. The edges of Q are formed by the intersection 
of the faces of C\ with those of C 2 . Therefore we need to also decompose the faces of 
C\. This is carried consistent with another, independent, percolation problem on an fee 
sublattice C' 2 of C 2 , and its dual C 2 . To each double decomposition of the faces of C\ into 
closed non-intersecting surfaces, and similarly of the faces of C 2 , corresponds a unique 
decomposition of Q into closed loops. 
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Figure 12: A decomposition of Q on the 3d L-lattice, showing how the loops (blue) reflect 
off the open bonds of percolation on the fee lattices C[ (solid magenta lines) and C' 2 (solid 
turquoise lines). 

This is the classical loop model which we seek. Unlike the 2d case, it corresponds to 
two independent percolation problems, one on an fee lattice of Ci, the other on an fee 
sublattice of C 2 . Each closed loop of Q is formed by the intersection of a closed surface 
made up of faces of C\ and a closed surface made up of faces of Ci. (It is of course possible 
that such pairs of closed surfaces intersect in more than one loop of Q, or not at all.) It 
alternately 'reflects' off an open edge of percolation on C[ (or an open dual edge of C"), 
then an open edge of C 2 or of C 2 , and so on (see Fig. IT21 . 

We now discuss the physics of this model and its implications for the class C network 
model. In principle we can assign different probabilities, p\ and P2, to the two independent 
percolation problems on C\ and Ci- The phase diagram is symmetric under the duality 
symmetry pj — > 1 — pj, so we can restrict attention to the quadrant p\ < |, p 2 < 2 - For 
either p\ or p 2 less than the bond percolation threshold p^ c ~ 0.12 [TH], since the clusters 
are finite so are the hulls on either C\ or C 2 (or both), and, since the loops on Q are 
formed by their intersection, these must be finite also, corresponding to a localized phase. 
Preliminary simulations of the model[20] for values of pj close to | appear to show that 
the loops are no longer finite, but, unlike the case of the diamond lattice [18J, neither are 
they simple random walks on large scales with fractal dimension 2, as would be expected 
of a sample exhibiting ohmic behavior. In fact, their fractal dimension appears to be close 
to 3, indicating that they are space filling. This may be a pathology of this model. It is 
known that the fractal dimension of cluster boundaries for p > p c is 3 - that is, a finite 
fraction of the infinite cluster, which has dimension 3 because it contains a finite fraction 
of the sites, is on its boundary. The loops for p > p c are formed by the intersection of 
random such df = 3 objects, so it is perhaps not surprising that they should also have 
d f = 3. 

For this reason it may be more useful to consider a 3d version of the Manhattan 
lattice. This is based on the same graph Q, but the edges are oriented so that along each 
line through the lattice they point in the same direction, and the lines alternate direction 
transversally. As in the 2d case, the probability of turning at a given node is p and of 



18 



going straight on 1 — p. (In principle we could again take different values of p at the two 
types of nodes corresponding to C\ and C 2 .) Analogously with the 2d case (see Sec. I4.2p . 
we can introduce an associated percolation problem on fee sublattices of C2 and C±, so 
that if, in the decomposition of Q, the paths turn at a given node, they reflect off an open 
edge. This means that, as in 2d, they are constrained to lie on regions occupied by the 
dual clusters of each percolation problem. Thus, for p > 1 — p^ c w 0.88, the loops of Q 
are almost surely finite in length, corresponding to the existence of a localized phase of 
the corresponding class C network model. 

On the other hand, for small p we expect a finite fraction of the paths to escape to 
infinity on a infinite lattice, corresponding to extended states. In this case the particle 
following a path will, almost, all the time, go straight ahead, with only a small probability 
p <C 1 of turning. In this case, the whether the particle has traversed an even or odd 
number of edges should be unimportant, leading to an effective simple random walk with 
diffusion constant 0{p^ x ) on intermediate distance scales. On larger scales, the walk may 
revisit regions it has in the past, but, unlike the case of 2d, this is unlikely because the 
3d random walk is not recurrent. However, as appealing as this argument may be, it is 
not rigorous, and indeed its proof appears to be of the same order of difficulty as showing 
that the 'true' self-avoiding walk in 3d is asymptotically gaussian|21[ 122] . Thus, at this 
stage, a proof of the existence of extended states in this class of network models remains 
elusive. 

6 Summary and further remarks 

We have shown how the quantum — > classical mapping for class C network models helps 
to gain insight into the nature of Anderson localization in both two and three dimensions. 
These models have direct physical relevance in systems where time-reversal symmetry is 
broken but spin rotational symmetry is preserved. The classical models correspond to 
deterministic motion in a random medium with two-sided mirrors, or, equivalently, cer- 
tain kinds of history-dependent random walks. On some lattices, including the important 
examples of the L and Manhattan lattices, the paths are associated with the hulls of per- 
colation clusters, and rigorous information can be inferred on whether the corresponding 
quantum model is in a localized or extended phase. 

However, there are a number of unresolved questions. Although the diamond lattice 
and 3d Manhattan lattice discussed in Sec. [5] are expected to exhibit an Anderson transi- 
tion, there is as yet no proof of the existence of an extended phase in which the motion is 
asymptotically diffusive, although this is strongly indicated on both numerical and other 
grounds. The relation of these models to other types of history-dependent random walks, 
such as the 'true' self-avoiding walk, in which the walk avoids regions it has visited in the 
past, is also unclear. Like Anderson localization, d = 2 is a critical dimension for the true 
self-avoiding walk. However in this case the RG flows |22[ |2"3"] are to free random walks for 
d > 2, and to a non-trivial stable fixed point for d < 2, while for Anderson localization 
we expect to find a non-trivial unstable fixed point for d > 2. This suggests that the 
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two problems are related by a change of sign of the interaction. However the analysis of 
Peliti and Obukhov[23] shows that for history-dependent random walks there are in fact 
three coupling constants which are potentially important near d = 2. An attempt to fit 
the walks on the 2d Manhattan lattice into this picture was made in Ref. [IS]. However 
this was not systematic and further work needs to be done in this direction. It should 
be noted, however, that a sigma-model analysis of the original class C quantum model 
does give the expected unstable fixed point for d > 2. [24]. A related question is that of 
the upper critical dimension for the transition in this model. The relation to interacting 
random walks suggests that this might be d = 4, as for ordinary polymers, but in this 
case the interactions are not simply repulsive, so this conclusion may not hold. 
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